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Abstract 

We numerically simulate the low-energy properties of interacting electrons in a random potential using the Hartree- 
Fock based exact diagonalization method. In particular, we investigate how the transport properties are influenced 
by the combined effects of disorder and correlations in the presence of the electron spin. To this end we calculate 
the participation number of many-particle states in Fock space, the return probability of single-particle excitations, 
and the Kubo-Greenwood conductance. It turns out that in the strongly localized regime interactions increase 
the conductance whereas for weak disorder interactions decrease the conductance. In contrast, single-particle 
excitations in general experience a localizing influence of the interactions. 



Key words: electronic transport, disorder, localization, metal-insulator transitions 



1. Introduction 

The interplay between Anderson localization 
and electron-electron interactions belongs to the 
most important unsolved problems of today's con- 
densed matter physics. The standard theoretical 
approach to this problem is based on perturba- 
tion theory in both disorder and interactions. In 
the metallic regime, the interactions lead to quan- 
tum corrections to the conductivity in addition 
to the weak- localization corrections [1] . A scaling 
theory [2,3] was developed on the basis of these 
perturbative results. It predicts that a phase tran- 
sition between an insulator and a normal metal, 
i.e., a disordered Fermi liquid, is possible only for 
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dimensions d > 2. In d = 2 the situation is incon- 
clusive. Under renormalization the disorder flows 
to zero but the interaction strength diverges, i.e., 
the perturbation theory becomes uncontrolled. 

The attention to these questions resurged in 
the mid nineties, mainly due to two surprising 
developments. First, Shepelyansky suggested [4] 
that in the localized regime two interacting par- 
ticles can form a pair whose localization length is 
much larger than that of a single particle. Second, 
measurements [5] on high-quality Si-MOSFETs 
revealed indications of a true metal-insulator 
transition (MIT) in two dimensions. Since these 
experiments are carried out at low electron den- 
sity where the Coulomb interaction is particularly 
strong compared to the Fermi energy, interaction 
effects arc a likely reason for this phenomenon. A 
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complete understanding has, however, not yet been 
obtained. In a number of theories new quantum 
states of matter have been postulated, either non- 
Fermi liquid metals or unusual superconductors 
[6-8]. However, other explanations based on more 
conventional physics like temperature-dependent 
disorder [9] or impurity screening [10], imply that 
the seeming MIT is a transient phenomenon, and 
the true ground state is always an insulator. 

We have recently developed [11,12] an efficient 
method to investigate disordered interacting elec- 
trons numerically, the Hartree-Fock (HF) based 
exact diagonalization (HFD). It is related to the 
quantum-chemical configuration interaction ap- 
proach. We have used this method to study the 
influence of interactions on the conductance of 
spinless fermions in one, two, and three dimen- 
sions. We found a delocalizing tendency of the 
interactions for strong disorder but a localizing 
one for weak disorder. In this work we extend our 
HFD investigations to electrons with spin. In par- 
ticular, we calculate the participation number of 
the many-particle eigenstates in Fock space, the 
localization properties of single-particle excita- 
tions, and the Kubo-Greenwood conductance. 

2. Model and method 

The generic model for spinless interacting dis- 
ordered electrons is the quantum Coulomb glass 
[13,14]. In this paper we use a generalization of the 
quantum Coulomb glass to electrons with spin. It 
is defined on a regular square lattice with M = L 2 
sites occupied by TV = N-\ + N± — 2KM electrons 
(0 < K < 1). To ensure charge neutrality each lat- 
tice site carries a compensating positive charge of 
2Ke. The Hamiltonian is given by 

H=-t^2 ( c L c i<? + h - c -) + V ini ° ( 1 ) 

{ij),a i,a 



where cJ CT , c, CT , and n ia are the creation, annihila- 
tion and occupation number operators at site i and 
spin a. (ij) denotes all pairs of nearest-neighbor 
sites, t is the strength of the hopping term, i.e., 
the kinetic energy. We parametrize the interac- 
tion Uij — e /rij by its value U between nearest- 
neighbor sites and include a Hubbard interaction 
Uh- The random potential values ipi are chosen in- 
dependently from a box distribution of width 2Wq 
and zero mean. The boundary conditions are peri- 
odic and the Coulomb interaction is treated in the 
minimum image convention (which implies a cut- 
off at a distance of L/2). 

A numerically exact solution of a quantum 
many-particle system requires the diagonalization 
of a matrix whose dimension increases exponen- 
tially with system size. This severely limits the 
possible sizes. To overcome this problem we have 
developed the HFD method. The basic idea is to 
work in a truncated Hilbert space consisting of the 
HF ground state and the low-lying excited Slater 
states. For each disorder configuration three steps 
are performed: (i) find the HF solution of the prob- 
lem, (ii) determine the B Slater states with the 
lowest energies, and (iii) calculate and diagonalize 
the Hamilton matrix in the subspace spanned by 
these states. The number B of basis states deter- 
mines the quality of the approximation, reasonable 
values have to be found empirically. 



3. Results 

To characterize, to what extent the interactions 
introduce non-trivial correlations (beyond the HF 
level) into the system we calculate the participa- 
tion number of the many-particle ground state 
with respect to the HF Slater determinants. The 
results for a system of 4 x 4 sites at half filling are 
shown in Fig. 1. For U = the Fock space partici- 
pation number is 1 because the HF approximation 
is (trivially) exact for non-interacting electrons. 
Increasing the interaction mixes the single-particle 
states, and the inverse Fock space participation 
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Fig. 1. Inverse Fock space participation number fj^, k as a 
function of interaction U for a system of 4 X 4 sites occupied 
by 8 spin-up and 8 spin-down electrons. The disorder is 
W = 1, the Hubbard energy U H = U, and the HFD 
basis size B = 500. The data represent averages over 1000 
disorder configurations. 

number is reduced. Nontrivial interaction effects, 
if any, become manifest mostly in this region. For 
even larger interactions the eigenstates become 
again more localized in Fock space since the HF 
approximation becomes exact again for U — > oo. 

Analogous information can be extracted from 
the single-particle survival probability 

= if E 1™ 5 - G *> + i6 ) G > ( 2 ) 

ij 

which averages the sum of the transmission prob- 
abilities of a single-particle excitation to all sites. 
Here G^-' A (u) are the retarded and advanced 
single-particle Greens functions. The survival 
probability at the Fermi energy, shown in Fig. 2, is 
related to the square of the quasiparticle weight. 

We now turn to the single-particle localization 
properties which we characterize by the return 
probability of the single-particle excitations, 

= E J™ 5 - G ^ + i5 ) G ^ i5 )- ( 3 ) 

i 

However, in this quantity localization and decay of 
the quasiparticles are entangled. To extract the lo- 
calization properties we normalize R p by the sur- 
vival probability Z. The data in Fig. 3 show that, 
in general, interactions tend to localize the single- 
particle excitations. This is in agreement with ear- 
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Fig. 2. Single-particle survival probability at the Fermi 
energy, Z(0); parameters as in Fig. 1. 
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Fig. 3. Normalized single-particle return probability at the 
Fermi energy, R p (0)/Z(0); parameters as in Fig. 1. 

lier HF results [14] which showed that the Coulomb 
gap in the single-particle density of states is respon- 
sible for the trend towards localization. In Fig. 3 
a tiny derealization seems to occur at low inter- 
action strength and small kinetic energy. However, 
the changes are within the statistical errors of our 
study. 

In a real transport experiment the most accessi- 
ble observable is the conductance. Theoretically, it 
can be obtained from linear-response theory. It is 
essentially determined by the current-current cor- 
relation function of the ground state. The real (dis- 
sipative) part of the conductance (in units of e 2 /h) 
at frequency uj is given by the Kubo-Greenwood 
formula [15], 
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Fig. 4. d.c. conductance G(0) for the same system as in 
Fig. 1. The inhomogeneous broadening is 7 = 0.0625. 

Re G -( w ) = HZE! V \(0\j*\v)\ 2 6(w + E - E U ).{A) 
(jj — 

j x is the x component of the current operator and 
\v) is an eigenstate of the Hamiltonian. Eq. (4) 
describes an isolated system while in a real d.c. 
transport experiment the sample is connected to 
contacts and leads. This results in a finite life time 
t of the eigenstates leading to an inhomogeneous 
broadening 7 = r _1 of the 5 functions in (4) [16]. 
To suppress the discreteness of the spectrum of a 
finite system, 7 should be at least of the order of 
the single-particle level spacing. Here, this requires 
comparatively large 7 > 0.05. We tested different 
7 and found that the conductance values depend 
on 7 but the qualitative results do not. 

The typical d.c. conductance is shown in Fig. 4. 
Since the logarithm of the conductance rather than 
the conductance itself is a self-averaging quantity 
in a disordered system we calculate the typical con- 
ductance by averaging the logarithms of the con- 
ductance values of 1000 disorder configurations. 
We find that, as in the spinless case, there are 
two qualitatively different regimes. For large ki- 
netic energy t the interactions always reduce the 
d.c. conductance, while for small kinetic energy, 
i.e. in the localized regime, small and moderate 
interactions significantly enhance the d.c. conduc- 
tance. For larger interaction strength the conduc- 
tance drops, indicating the crossover to a Wigner 
crystal or Wigner glass. 



A comparison with analogous results for spinless 
fermions [11] shows that the delocalizing influence 
of the interactions is significantly stronger for elec- 
trons with spin. 
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